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Orbiter Landing Loads Math Model Description and Correlation With 

ALT Flight Data 

David A. Hamilton , John A. Sch liesing, and George A. Zupp, Jr. 


SUMMARY 

In recent years, there has been a trend to use ra- 
tional analysis in defining the touchdown dynamics 
of aircraft. This trend has precipitated efforts to cor- 
relate these analyses with measured landing dynamic 
data. Along this line, a mathematical model of the 
Space Shuttle Orbiter landing dynamics was 
developed and verified using measured landing loads 
data from the Shuttle Approach and Landing Test 
(ALT) Program. The correlation between analysis 
and test demonstrates the ability to accurately pre- 
dict landing dynamics and the resulting landing loads 
on a flight vehicle. The key elements in the accuracy 
of correlation in the mathematical model are the 
landing gear strut stroking, fore-and-aft and lateral 
strut bending, and tire "spin-up” characteristics. Also 
modeled are the aerodynamic forces and moments 
on the Orbiter, including control surface effects. 

Detailed correlations in the time domain between 
analysis and test of the landing gear "spin-up” and 
"spring-back” loads are presented, with peak loads 
agreeing well in both magnitude and frequency. 
Other landing gear parameters such as strut stroke 
and wheel rotational velocity are also correlated. 
Conclusions from these studies are included. 


INTRODUCTION 

The advent of the Space Shuttle Program has add- 
ed a new dimension to the concept of reusable 
spacecraft. The Space Shuttle launch vehicle is com- 
posed of four principal structural elements: two solid 
rocket boosters, an Orbiter, and an external tank that 
provides liquid hydrogen and oxygen to the Orbiter 
main engines. The solid rocket boosters and the Or- 
biter are reusable. The primary sequencing points in 
the Shuttle trajectory are (1) lift-off, (2) solid rocket 
booster staging at approximately 120 seconds into 


the mission, and (3) end-burn and external tank stag- 
ing at approximately 500 seconds into the mission. 
At this point, the Orbiter continues on into Earth or- 
bit carrying payloads of as much as 29 500 kilograms 
(65 000 pounds). Once the payloads have been 
deployed, recovered, etc., and the orbital activity has 
been concluded, the Orbiter will enter the Earth’s at- 
mosphere and land. 

The Orbiter is a delta-wing spacecraft designed to 
enter the atmosphere at orbital speeds and maneuver 
to a designated landing site. The Orbiter lands un- 
powered on a conventional runway. The landing 
system is a typical tricycle configuration with a nose 
gear and two main gears. The landing system is 
designed to handle approach speeds above 100 m/s 
(200 knots) with a sink speed of as much as 2.9 m/s 
(9.6 ft/s). 

The magnitude of the initial landing loads dictates 
the size of a significant portion of the Orbiter struc- 
tures; therefore, precise knowledge of the landing 
loads can result in a significant structural weight sav- 
ing. In an effort to minimize the Orbiter structural 
weight, a detailed mathematical model of the Orbiter 
touchdown dynamics was developed and programed 
on a digital computer. 

In the detailed simulation of the Orbiter 
touchdown dynamics, the Space Shuttle Orbiter is 
idealized as a pseudorigid structure to which a land- 
ing gear system is attached. The landing gear system 
is made of telescoping oleo-pneumatic struts and 
pneumatic tires. The aerodynamic forces acting on 
the Orbiter during landing are approximated by "six- 
degree-of-freedom” nonlinear aerodynamic equa- 
tions with aerodynamic control coming from the 
deflection of the Orbiter’s rudder, body flaps, 
elevons, and ailerons. 

This report includes a definition of the pertinent 
coordinate systems and coordinate transformations. 
The equations of motion governing both a rigid vehi- 
cle and landing gear dynamics are presented. Con- 


siderable detail is provided in defining the kinematic 
relationships used to model the articulating landing 
gear geometry. 

The purpose of this report is twofold: (1) to docu- 
ment the analytical procedures used in developing a 
mathematical model of the Space Shuttle Orbiter 
touchdown dynamics and the quality of comparisons 
between the predicted results and the flight data, and 
(2) to make these procedures available to others who 


may be confronted with the task of mathematically 
modeling the touchdown dynamics of an oleo- 
pneumatic-type landing gear system. 

In compliance with NASA’s publication policy, 
the original units of measure have been converted to 
the equivalent value in the Systeme International 
d’Unites (SI). As an aid to the reader, the SI units are 
written first and the original units are written 
parenthetically thereafter. 


SYMBOLS 


^piston,- 

piston area of /'th strut, cm 2 (in 2 ) 

Fbx, 

bending force of /th strut in X B 
direction, N (lb) 

Qx, 

fore-and-aft damping coefficient of 
/th gear strut, N-s/m (lb-s/ft) 

F Byj 

bending force of /th strut in Y B 
direction, N (lb) 

C Byj 

lateral damping coefficient of /th 
gear strut, N-s/m (lb-s/ft) 

^bend; 

bending force in /'th strut, N (lb) 

Qlamp, 

damping coefficient of /'th strut 

F D, 

drag force on /'th gear perpendicular 

(nonlinear), N-s 2 /m 2 (lb-s 2 /ft 2 ) 

to strut axis and in X&-Z B plane, N 
(lb) 



[C)j 

matrix of y'th measurement calibra- 
tion constants 

Fdampf 

damping force in /‘th strut, N (lb) 

F>A\ r D B\j 

trunnion drag forces on /th gear, N 

F strut/ 

total strut force acting in stroking 

(lb) 


direction of /'th strut, N (lb) 

EI X , 

average fore-and-aft bending stiff- 

F Vi 

vertical ground reaction force of /'th 

ness of /th strut, N-cm 2 (lb-in 2 ) 


gear, N (lb) 

El yj 

average lateral bending stiffness of 

fx a ,fy a ,fz a 

aerodynamic force vector compo- 

/th strut, N-cm 2 (lb-in 2 ) 


nents in the body coordinate 
system, N (lb) 

Fa 

aerodynamic force vector at 
reference point, N (lb) 

hi 

side force on /'th gear perpendicular 
to strut axis and in the Y B direction, 
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aerodynamic force vector at c.g., N 


N (lb) 


(lb) 

hi 

polar moment of inertia of /'th 

F 

1 air/ 

force from compression of air 
spring of /'th strut, N (lb) 

wheel, kg-m 2 (slug-ft 2 ) 




hxJyyJzz 

mass moments of inertia of the Or- 

F *i 

lower bearing force on /'th gear, N 

biter, kg-m 2 (slug-ft 2 ) 

(lb) 

hzJyzJzz 

mass products of inertia of the Or- 


F B2 i 

upper bearing force on /'th gear, N 

biter, kg-m 2 (slug-ft 2 ) 

(lb) 

K bi 

axial stiffness of /'th gear drag brace, 


Fb, 

bearing friction force from stroking 

N/m (lb/ft) 


of /'th strut, N (lb) 


2 


gas constant for /th strut 

n 

rolling radius of /th tire, m (ft) 

damping coefficient of /th tire, 

s 

reference area, m 2 (ft 2 ) 

N-s/m (lb-s/ft) 

t a 

aerodynamic moment vector at 

spring rate of /th tire, N/m (lb/ft) 


reference point, J (ft-lb) 

length between /th gear bearing 

^c.g. 

aerodynamic moment vector at c.g.. 

points when gear is fully extended, 

J (ft-lb) 

m (ft) 

[t b A 

Euler angle transformation matrix 

strut length from axle to upper bear- 


from body coordinate system to in- 

ing on /th gear, m (ft) 


ertial coordinate system 


[Teb] 

transformation matrix from Euler 
angular rates to body angular rates 

fully extended strut length of /th 
gear measured from trunnion attach 

TX a JY a JZ a 

aerodynamic moment vector com- 

point to gear axle, m (ft) 


ponents at Orbiter c.g. in the body 
coordinate system, J (ft-lb) 

distance between upper strut trun- 
nion and drag brace upper trunnion 

TxJyJz 

components of total torque on Or- 

measured parallel to strut axis for 

biter, J (ft-lb) 

/th gear, m (ft) 

t 

time, s 

Mach number 

'/ 

time when the wheel speed is equal 

mass of the Orbiter, kg (slugs) 


to y a p s 

drag brace load for /th gear, N (lb) 

to 

initial time of touchdown, s 

load vector from strain gage reading 

V 

magnitude of the Orbiter velocity, 

for /th measurement, N (lb) 


m/s (ft/s) 

initial pressure in /th strut, N/m 2 

Ya, 

translational velocity of /th wheel in 

(lb/ft 2 ) 


X l direction, m/s (ft/s) 

aerodynamic pressure, N/m 2 (lb/ft 2 ) 

Y A 1 Y B \ i 

trunnion vertical forces on /th gear, 
N (lb) 

vector location of Orbiter c.g., m (ft) 

Ya 

local speed of sound, m/s (ft/s) 

vector location of axle of /th gear 
from the Orbiter c.g., m (ft) 

^ head; 

head volume of /th strut when fully 
stroked, m 3 (in 3 ) 

undeflected radius of /th gear tire, m 
(ft) 

V j. 

s i 

slip velocity of /th tire, m/s (ft/s) 

vector location of aerodynamic 
reference point, m (ft) 

Y tire,- 

velocity vector of /th tire, m/s (ft/s) 

vector location of upper trunnion of 
/th gear from the Orbiter c.g., m (ft) 

Y x ,Yy,V z 

Orbiter velocity components in 
body coordinate system, m/s (ft/s) 


fV x ,W y ,tV z 

ground-wind velocity components 

orientation vector of /th tire, m (ft) 

in body coordinate system, m/s 

(ft/s) 
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X, Y,Z Orbiter inertial accelerations, m/s 2 

(ft/s 2 ) 

A orthogonal axes of body coordinate 
system 

A 7,K/,Z/ orthogonal axes of inertial coordi- 

nate system 

a angle of attack, deg 

/3 sideslip angle, deg 

8 a aileron position, deg 

8 bp body-flap position, deg 

8 e elevon position, deg 

8 r rudder position, deg 

§5/ stroke of /th gear strut, m (ft) 

h s . stroking velocity of /th strut, m/s 

(ft/s) 

8$ n nose gear strut stroke, m (ft) 

8s r right main gear strut stroke, m (ft) 

8,. deflection of /th tire in Z t direction, 

m (ft) 

8 v . fore-and-aft deflection of /th strut in 

•*/ 

^direction, m (ft) 


0/ 

alinement angle between /th tire 
orientation and velocity vector, deg 

ft/ 

angular velocity of /th wheel about 
T^axis, rad/s 

(i) X ,U)y,(X) Z 

angular velocity components of Or- 
biter, rad/s 

Operators: 


o 

first derivative with respect to time 

(••) 

second derivative with respect to 
time 

(-) 

vector quantity 

(') 

body prime coordinate 

("> 

body double-prime coordinate 

("') 

body triple-prime coordinate 

n 

matrix 

n- 

matrix inverse 

u T 

matrix or column vector transpose 

{ ) 

set of vectors 

o 

unit vector 



lateral deflection of /th strut in Y B 
direction, m (ft) 

My 

strain gage vector for ./'th measure- 
ment 

0 xfiyfi Z 

Euler angles from inertial to body 
coordinate system, deg 


bearing friction coefficient of /th 
strut 

M/ 

drag force friction coefficient of /th 
tire 

Vsi 

side force friction coefficient of /th 
tire 

p 

local air density, kg/m 3 (slugs/ft 3 ) 


MATHEMATICAL IDEALIZATION 

Coordinate Systems and Transformations 

Two coordinate systems are used in the formula- 
tion of the mathematical simulation of the Shuttle 
vehicle touchdown dynamics. Each coordinate 
system is an orthogonal right-hand set. 

Inertial coordinate system (X t , Y t , and Z { axes ). — 
The inertial coordinate system is defined such that 
the X r Y/ plane forms the landing surface and is 
assumed to be universally fixed (fig. 1). The Z 7 axis is 
perpendicular to and directed up from the X r Y I 
plane. The gravity vector is taken to be parallel to the 
inertial Z 7 axis. The equations of motion of the vehi- 
cle and landing gear are expressed in the inertial 
coordinate system. 
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Body coordinate system (X B , Y B , and Z B axes ). — 
The body coordinate system is fixed in the Shuttle 
vehicle so that the origin coincides with the vehicle’s 
center of mass. The body X B axis is directed along 
the longitudinal geometric centerline of the vehicle 
as illustrated in figure 1. 

Vehicle geometry, such as landing gear trunnion 
points, wheel locations, etc., is described by vectors 
expressed in the body system. The body coordinate 
system is oriented with respect to the inertial coordi- 
nate system by three Euler angles. These angles are 
produced by three successive rotations about the 
three body axes in a prescribed order. Associated 
with a given rotation order is the Euler angle 
transformation, which transforms vectors in the 

body system to vectors in the inertial system. This 
transformation is expressed in matrix form as 


(*/] 


[ X B 

Y I 

\Jbi] 

y b 

\ Z I J 


Z B 


where the elements X J} T 7 ,Z/and X B , Y B ,Z B are the 
inertial and body coordinates, respectively. 

The Euler angle rotation sequence that will be 
used in the development of the Euler angle transfor- 
mation will be successive rotations about the body 
X B axis (9 X roll), Y-axis (0 y pitch), and Z^axis (0 Z 
yaw), respectively. The detailed derivation of the 
Euler transformation and Euler rate transformation 
is presented in the appendix. 


Equations of Motion 

The landing configuration of the Shuttle vehicle is 
idealized as a rigid body of constant mass. Attached 
to the vehicle are three landing gear assemblies, a 



FIGURE 1. — Coordinate system definition. 


nose gear and two main gears (fig. 2). The mass of 
each landing gear assembly is assumed to be a point 
mass concentrated at the wheel axle of each assem- 
bly. Instead of using the actual mass of the landing 
gear assembly as the wheel mass, an “effective wheel 
mass” is computed so that the effective wheel mass 
located at the wheel axle approximates the inertial 
characteristics of the landing assembly. Similar 
assumptions have been made in studies reported in 
references 1, 2, and 3. 

The rigid vehicle has six degrees of freedom, three 
in translation and three in rotation. Each landing 
gear assembly has three degrees of freedom in 
translation and one degree of freedom in rotation. 
The three translational equations of motion of the 
vehicle center of mass are 


X 

_ j_ 

ZFX 

Y 

ZFY 

Z 

M o 

ZFZ 


which is a mathematical expression of Newton’s sec- 
ond law. The three rotational degrees of freedom are 
described by the three Euler angles (0 X ,0 V ,9 Z ) and are 
related by the three rotational equations of motion 
(equations (44), (45), and (46) in the appendix). 

The rotational equations of motion are functions 
of the angular rate and angular acceleration vectors 
and the vehicle inertia. The angular acceleration vec- 
tor is determined by solving the Euler moment equa- 
tions or 


(T \ 
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Using an appropriate integration scheme, the solu- 
tion of the vehicle equation of motion will yield the 
time histories of the characteristic parameters. 
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(a) 

FIGURE 2. — Landing gear chassis assembly. 


Landing Gear Forces 

The Orbiter landing gear configuration is a tricycle 
type with primary braking on the main wheels. Dur- 
ing landing impact, the vertical motion of the Orbiter 
is arrested through the compression of the oleo- 
pneumatic struts and tires of the landing gear system. 

The landing gear forces that are critical to the Or- 
biter structural design are those associated with the 
landing gear wheel “spin-up” and “spring-back” 
phenomenon that occurs during gear impact. As the 
aircraft approaches the landing site, the rotation 
speed fit, of the gear wheels is approximately zero. 
As the wheels make contact with the landing surface, 
the friction between the tires and the surface causes 
the wheels to rotate. The frictional forces increase 
the rotational speed of the wheel until it is syn- 
chronous with the translational velocity of the gear 
axle V / 4 r 

Because of the elasticity of the gear, the frictional 
forces will cause the landing gear strut to deflect. At 
some point in the aft deflection of the strut, the land- 
ing gear trunnion forces will be at a maximum. This 
point is defined as the maximum spin-up load. The 
subsequent rebound of the strut will cause a reversal 



Main landing gear, (b) Nose landing gear. 


in loading in the gear upper trunnions. This reversal 
in loading is defined as the spring-back load. Usually 
the spin-up and spring-back loads are at a maximum 
during the first bending oscillatory cycle. 

The landing gear forces during initial impact are 
governed by the tire compression modulus. During 
this phase, the tire load increases with deflection un- 
til the force component along the strut axis is greater 
than the preload in the strut, at which time the gear 
will begin to stroke. Prior to this, the oleo-pneumatic 
strut will remain fully extended. The magnitude of 
the strut preload is set by the initial internal pressure 
of the strut air spring. As the tire load exceeds the 
strut preload, the strut will start to stroke. The strut 
stroking forces are due primarily to the nonadiabatic 
compression of the air spring, the damping fluid 
moving through an orifice, and the bearing friction. 

The landing gear bending characteristics are ap- 
proximated by a massless elastic beam with damp- 
ing. The elastic beam connects the gear axle to the 
Orbiter gear trunnions. The inertia of the gear is ap- 
proximated by a mass located at the gear axle; this 
mass is defined as the unsprung mass. The idealiza- 
tion of the gear is pictorially presented in figure 3. 
The unsprung mass is given four degrees of freedom, 
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FIGURE 3. — Landing gear idealization. 


three in translation and one in rotation. The rota- 
tional degree of freedom is associated with the spin- 
up phenomenon of the wheel assembly_during the 
initial touchdown phase. The vectors R tj and R gj 
define the positions of the upper trunnion and gear 
wheel axle, respectively. 

The forces acting on the unsprung mass can be 
divided into ground reaction forces and strut forces. 
The orthogonal components of the ground reaction 
forces are (1) vertical force F Vj , (2) drag force Fd p 
and (3) side force F^.. The side and drag forces are 
assumed to be in the plane of the landing surface. 
Similarly, the strut forces are due to the strut strok- 
ing dynamics and strut bending. The strut bending 
forces are orthogonal to the longitudinal axis of the 
strut and act at the wheel axle. 

Tire/ground forces . — The tire vertical ground forces 
are a function of the tire deflection 8,. and deflection 
rate 8 r/ and are computed by 


\ = *tir ef'i + K « T S 


( 4 ) 


where A^ re/ is the nonlinear elastic modulus of the 
tire and K^f j \s the nonlinear damping coefficient of 
the tire. During wheel spin-up, the friction force act- 
ing on the tire is a function of the vertical load F v/ and 
the slip velocity V Sj . From experimental data, the 
friction coefficient /jl 1 -between the tire and the land- 
ing surface is developed as a function of the slip 
velocity. Generally, the friction coefficient is at a 
maximum when the wheel rotational velocity is ap- 
proximately 80 to 90 percent of the axle translational 
velocity. The drag force Fq. on the tire is computed 
by 




( 5 ) 


The slip velocity is given by the solution of the in- 
tegral equation 
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where I Pj is the polar moment of inertia of the land- 
ing gear wheels, r, is the rolling radius of the fth tire, 
t Q is the initial time of touchdown, and ^ is the time 
when the wheel speed is equal to V Aj . 

It is assumed that the tire side force F$. is a func- 
tion of the angular alinement between the tire orien- 
tation and the tire translational velocity and the ver- 
tical load. (It should be noted that this assumption 
for tire side forces was adequate for the approach and 
landing test (ALT) correlation. This model of tire 
side force is not representative when the tire aline- 
ment angle and vertical load are “large.” When very 
accurate tire side force models are required, models 
similar to those used in shimmy analysis should be 
used (ref. 4).) 

If the orientation of the tire is denoted by the vec- 
tor /? t j re/ and the tire translational velocity vector, in 
the plane of the landing surface, is denoted by 
then the alinement angle is computed by 


D, = cos 


-l 



-y 

R . 

• 

— y 

v. 

“ 


tirey 


tirey 



— y 


— y 



R . 

• 

v . 


__ 

tirey 


tirey 
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(7) 


ternal forces are produced by the nonadiabatic com- 
pression of the air spring, by the bearing friction, and 
by the damping fluid being forced through a metered 
orifice. 

The magnitude of the damping force is governed 
by the area of the orifice. The orifice area is con- 
trolled by a metering pin the diameter of which is a 
function of the strut stroke position. The design of 
the cross-sectional area along the length of the meter- 
ing pin is governed by the energy absorption require- 
ments of the landing gear. In the case of the Orbiter, 
the maximum design sink rate is 2.9 m/s (9.6 ft/s). 
At the initial stage of strut stroking, the stroke rate is 
of the same order as the sink speed. During this 
period, the air spring forces are low and the damping 
forces are high. As the strut continues to stroke, the 
air spring forces will increase and the damping forces 
will decrease such that the total force (air spring and 
damping) is approximately constant over the strut 
stroke distance. This strut load-stroke design will 
maximize the landing gear energy absorption 
capability for a given maximum allowable axial strut 
force. 

The strut damping forces ^damp/C 311 be computed 
by 


The tire force, acting normal to the tire orientation 
vector /? t i re/ , normally increases when the magnitude 
of the alinement angle increases to the point of tire 
skidding. 

The side force is also dependent on the magnitude 
of the tire vertical load F Vj . From experimental data, 
the side force friction coefficient fjL Sj can be expressed 
as a linear function of the tire alinement angle 0 y. 
Under this assumption, the side force F^. is com- 
puted by 


^damp. ^damp- ^S- ( • 

* 1 1 6 


where 85 . is the strut stroking velocity and Q amp/ is 
the nonlinear damping coefficient, which is a func- 
tion of the strut stroke. The strut air spring force F ^ . 
can be computed by 


F. - n F„ 

0. ^s. v. 

M l l 


( 8 ) 


The coefficient /jl s . is normally developed from ex- 
perimental data and is presented as a function of the 
alinement angle </>,-. The side friction coefficient n Sj is 
approximately linear with the alinement angle to ap- 
proximately 10°. Alinement angles greater than 10° 
produce quasi-skidding conditions in the tires, result- 
ing in a reduction in tire side force. 

Strut forces . — The oleo-pneumatic strut attenuates 
the landing forces through the mechanism of strut 
compression. During the compression, the strut in- 


F -PA 
airy o. piston^. 


'~ K t 


heady 


heady 


A ~ 8 S. 

pistoriy 1 


K i 


( 10 ) 


where P 0j is the initial static internal strut pressure, 
^piston, is the strut piston cross-sectional area, K head . 
is the head volume of the strut when fully stroked, 
and Kj is the polytropic gas constant. 

The strut bending forces are reacted internally by 
two bearings: the lower bearing, which is fixed to the 
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upper section of the strut, and the upper bearing, 
which is fixed to the stroking segment of the strut. 
The bearing reaction forces F B \ j and F B2/ are illus- 
trated in figure 4. The bearing reaction forces are 
computed by 


The net axial strut force F strut/ is the summation 
of the components or 


struts 


F Br i + F aii i + ^damp. 


(14) 



(ID 


and 





* 

1 


( 12 ) 


where L, is the strut length from the axle to the upper 
bearing and L^ j is the length between the two bear- 
ings when the gear is fully extended. 

The bearing friction force F Brj present when the 
gear is stroking is given by 


% = \ 


If l + 1 f i 
\ bi^ y B 2J 

i 



hi 


(13) 


where /jl Bi is the bearing friction coefficient at the up- 
per and lower strut bearings. It should be noted that 
these bearing friction forces are present only in a 
dynamic condition; special care should be taken 
when programing these equations for the digital 
computer to ensure that these forces are not com- 
puted while in a static condition. 



The strut bending force F bend . acts normal to the 
strut axial force and is located at the wheel axle. The 
bending forces are divided into side-to-side bending 
and fore-and-aft bending. The magnitude of the strut 
bending forces can be related to the strut bending 
stiffness, bending deflection, and deflection rate. The 
bending deflection from side to side and the deflec- 
tion rate will be denoted by 8 y . and 5 V/ , respectively. 
Similarly, the fore-and-aft deflection and the deflec- 
tion rate are 8 X . and 8 X/ ., respectively. 

The fore-and-aft bending motion of the gear can 
be described by idealizing the gear, in the x direction, 
as a beam supported by a pinned joint at the upper 
trunnion point and a pinned joint at the drag brace 
attachment. An effective bending stiffness EI Xj is 
assumed. It is also assumed that the drag brace stiff- 
ness can be approximated by a linear spring of stiff- 
ness Kb r Under these assumptions, the fore-and-aft 
bending force Fg Xj is 


F Bx. C BxPx + 


/ 2 r 3 

N i l i 

( l n* - 2L n] + l n) 

K b. 3£ V 

_ 

_ 


(15) 


where 



(16) 


L 0j is the initial strut length, L\ j is the distance be- 
tween the strut upper trunnions and the drag brace 
attachment to the strut, C^ X/ is the fore-and-aft 
damping coefficient, and EI Xj is the average fore-and- 
aft bending stiffness. 

The side-to-side bending motion of the gear can be 
described by idealizing the gear, in the y direction, as 
a cantilever beam with an effective bending stiffness 
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EI yj . Under these assumptions, the bending force in where 

the y direction is 




(17) 


V‘ Vci- “y 2 ♦ OJ - "M 2 + K - “i ) 2 

( 21 ) 


The Mach number M is defined by 


where Cg yj is the side-to-side damping coefficient. 
The net bending force Fbend used in the calculation of 
the bearing forces is 



( 22 ) 


n«n d , “ + V, 2 < 18) 

With the forces acting on the gear unsprung rr.ss 
defined by equations (4), (5), (8), (14), (15), and 
(17), the accelerations of the unsprung mass can be 
computed. Using an appropriate numerical integra- 
tion scheme, the time history of the position vector 
R g can be determined. The landing gear forces acting 
on the Orbiter are assumed to be opposite and equal 
to those strut forces acting on the unsprung mass. 


Aerodynamic Forces 

Aerodynamic forces and moments acting on the 
Orbiter during landing are computed as a function of 
angle of attack, sideslip angle, Mach number, 
dynamic pressure, reference area, control surface 
positions, and aerodynamic coefficients. 

The angle of attack a (ref. 5) is computed from 
the Orbiter velocity vector and the local winds; i.e., 



where Kis the magnitude of the Orbiter velocity and 
V a is the local speed of sound. 

The dynamic pressure ~q is computed by 

q = \ P V 2 (23) 


where p is the local air density. 

The Orbiter autoland control system is modeled in 
this analysis according to its description in reference 
6. The control system determines control surface 
deflections as a function of Orbiter attitudes, angular 
rates, dynamic pressure, runway position, load fac- 
tors, horizontal velocity, and gear deflection. The 
surfaces include the elevons, ailerons, rudder, speed 
brakes, and body flaps. 

Aerodynamic data for the Orbiter landing con- 
figuration are described in reference 7. These non- 
linear data are used in these analyses in tabular form 
as a function of the aforementioned parameters. 

Thus, the six-degree-of-freedom aerodynamic 
forces and moments (FX A , FY A , FZ A , TX A , TY A , and 
TZ a ) are computed by equations of the form 

FX A = CxqS (24) 


where V z and V x are components of the Orbiter 
airspeed in the Orbiter coordinate system and W x 
and W z are the crosswind components in the same 
system. Similarly, the sideslip angle p is computed by 

0 = sin - 1 (-^ — p — (20) 


where S is the reference area and the aerodynamic 
coefficient is functionally of the form 

Cx = f(«,P,M,5 e ,b a ,8 r 5 SB ,5 BF ) (25) 

and the other forces and coefficients are of similar 
form. 
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The forces T A and the moments 7^ are transferred 
from the reference point to the center of gravity 
(c.g.) by 

?„ cg - f A ( 26 ) 

and 

\ t - t a * [K - *4 x %, < 27 > 

These are the aerodynamic forces and moments act- 
ing on the Orbiter that drive the Euler moment equa- 
tions and the Newton translational equations. 


GEAR INSTRUMENTATION 

The main and nose landing gears on Orbiter vehi- 
cle 101 were instrumented with strain gages on the 
lower drag link, the upper strut trunnion pins, and 
the strut/wheel torque link mechanism (fig. 5). The 
strain gages all had a frequency response of 0 to 50 
hertz. The four gages on the lower drag brace were 
mounted on an I-beam section of the brace and the 
output of the four gages was electrically averaged to 
obtain the resulting tensile or compressive strain. 
The strut trunnion pins were hollow cylinders with 
four gages mounted on the inner surface. The out- 
puts of the upper and lower gages on the pins were 
combined to obtain drag load and the outputs of the 
forward and aft gages were combined to obtain verti- 
cal load. Differential vertical loads between the in- 
board and outboard trunnions were used to obtain 



FIGURE 5. — Main and nose gear strain gage schematic. 


side loads. Each side of the main gear torque link 
mechanism was instrumented with a strain gage. The 
output of these two gages was combined to obtain 
strut torque. 

Each gear was also instrumented with a stroke 
measurement device and with wheel-speed sensors. 
All the above data were recorded on a wide-band 
recorder at a sample rate of 400 samples per second. 
This type of sample rate provided a good resolution 
of gear loading during the landing phase. 


LANDING GEAR CALIBRATION 

The main and nose landing gears were load cali- 
brated using a system calibration technique. 1 ’ 2 This 
calibration involved applying various combinations 
of vertical, drag, and side loads and torques to the 
gear at the axle as shown in figure 6 and then record- 
ing the strain gage outputs for each load condition. 
The load combinations were applied in incremental 
percentages of 20, 40, 60, 80, and 100 percent of max- 
imum load so that repeatability and linearity of the 
strain gage output could be established. 

The resulting set of load and output data was used 
in a “least-squares” solution to develop calibration 
constants for each gear. These calibration constants 
represent a linear relation between strain gage output 
and the corresponding gear loads. This relation is of 
the form 


{P} f = [C] f {e} f (28) 


where {e} 7 - are strain gage readings, [C] y are the 
calibration constants, and \P} ) are the corresponding 
loads. The constants [C] ; for the left main, right 
main, and nose gears are given in figure 7 and are 
taken from internal data. 1 ' 3 The strut was calibrated 


*A. L. Carter, “OV-101 Landing Gear Calibration Results,’' 
NASA Dryden Flight Research Center Internal Letter, July 1977. 

2 A. L. Carter, “OV-101 Landing Gear Post Flight Strain Gage 
Hookup and Sensitivity Checks," NASA Dryden Flight Research 
Center Internal Letter, Jan. 1978. 

3 Walter B. Horne, “Tire Friction Coefficient Estimates for 
Decelerating and Steering the Space Shuttle Vehicle During Land- 
ing on Conventional, Grooved, and Lakebed Surfaces," NASA 
Langley Research Center Internal Letter, Nov. 5, 1976. 
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Nose gear 


by stroking the strut in incremental amounts from 
fully extended to fully stroked and recording the 
electrical output of the linear potentiometer for each 
strut position. The wheel-speed sensors were calibra- 
ted in a laboratory but not on the gear assembly. 
However, their apparent accuracy was satisfactory. 


Right main gear 


’ db l 

= 43.272 

*9235 

P ° b n' 

19.65 

*9135 

V 

■ 76.621 

*9241 

V 

11.22 

*9141 

A 

• 76.832 

*9238 

V 

- 10.87 

*9140 

A 

= Not available 

° A V 

19.15 

*9139 

A 

= 40.294 

*9239 

° b1 n' 

: 19.43 

*9138 


LANDING LOADS AND CORRELATION 

Correlation with the ALT loads was done for both 
main and nose gear impact and for occurrence of 
maximum main gear tire loads. In general, the strain 
gage output from the right main gear was more con- 
sistent than that from the left main gear and thus was 
considered more representative of the actual loads. 

The horizontal velocity and sink rates at both 
main and nose gear impact for each ALT flight are 
shown in table I. The main gear loads for free flight 5 
(FF5) were the highest of all flights because of the 
high sink rate, whereas the nose gear loads were the 
highest for FF2. However, for nose gear correlation, 
FF5 loads were used because the FF5 flight repre- 
sented the actual Orbiter tail-cone-off configuration. 
The maximum main gear tire loads occurred on FF4 
because of the high horizontal velocity at nose gear 
impact, which produced a high downward aero- 
dynamic load on the Orbiter and thus a high tire load. 

As shown in the sections on tire/ground loads and 
strut loads, the forces acting on the gear strut at the 
axle can be calculated as a function of time. The cor- 
responding loads acting at the gear trunnions, as il- 
lustrated in figure 7, can be resolved from the loads 
at the axle using the gear geometry. The trunnion 
loads and the drag brace loads are used for correla- 
tion because they are directly proportional to the 
strain gage measurements. 


Left main gear 


v 

42.879 < 

9335 

V 

74.626 < 

9341 

V 

69.281 , 

e 9340 

V 

-142.29] 

1 *9339 

°B1 ■ 
B1 R 

' - m - 556 *9338 



FIGURE 7. — Main and nose gear calibration constants. 


The main gear impact loads on the right main 
lower drag brace and upper strut trunnion for FF5 
are shown in figure 8. The location of these loads is 
given in figure 7. Also shown with the measured 
loads are the predicted loads using the pre-ALT main 
gear math model. These three loads were used for 
correlation because they were largest in magnitude 
and all other ground reaction and trunnion loads can 
be calculated from them. The correlation between 
analysis and test was good for these loads. However, 
a refinement of strut bending damping and gear 
unsprung mass produced excellent correlation as 
shown in figure 9. The remaining differences be- 
tween analysis and test are attributed to inaccuracies 
in initial conditions, data measurement, and mathe- 
matical modeling. The right main strut stroke for the 


Trunnion f 



(a) (b> 

FIGURE 6.— Calibration test data, (a) Sketch showing applied 
loads, (b) Typical strain gage output (vertical load, inboard trun- 
nion, right main gear). 


Table I . — Horizontal and Vertical Velocities 
for ALT Landings 


Flight 

Main gear impact 

Nose gear impact 


Horizontal 

Sink rate. 

Horizontal 

Sink rate, 


velocity, 
m/s (knots) 

m/s (ft/s) 

velocity, 
mis (knots) 

m/s (ftls) 

1 

99 (192) 

~0.3 (-1.0) 

76 (148) 

1.3 (4.2) 

2 

96 (186) 

-.3 (-1.0) 

70 (137) 

2.1 (6.8) 

3 

99 (192) 

-.3 (-.9) 

77 (150) 

1.3 (4.3) 

4 

102 (199) 

.9 (3.0) 

84 (163) 

1.0 (3.3) 

5 

78 (152) 

a 2.0 (6.5) 

68 (132) 

2.0 (6.4) 


^Second impact. 
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FIGURE 8. — Right main landing gear loads and correlation 
using pre-ALT math model, (a) Drag brace, (b) Inboard trun- 
nion, vertical load, (c) Outboard trunnion, vertical load. 
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FIGURE 9. — Right main landing gear loads and correlation 
using updated math model, (a) Drag brace, (b) Inboard trunnion, 
vertical load, (c) Outboard trunnion, vertical load. 
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FF5 landing is shown in Figure 10. Correlation was 
excellent during the initial impact but there was 
about a 10-percent error at the time of maximum 
stroke. This error was also attributed to modeling 
and instrumentation inaccuracies. 

The wheel-speed spin-up transients for the right 
and left main gears for FF5 are shown in Figure 11. 
These data are shown in meters per second (knots) 
and represent the equivalent translational velocity 
that corresponds to the wheel angular velocity times 
the apparent rolling radius. The wheel speeds on 
each gear had initial values of approximately 30 m/s 
(60 knots) because of their partial spin-up on the pre- 
vious gear contact which occurred on FF5. Correla- 
tion was within 10 percent on magnitude and within 
20 percent on timing. It should be noted that the 
spin-up timing is very sensitive to the tire/ground 
friction coefficient. The friction curve used for the 



FIGURE 10.— Right main gear strut stroke for main gear impact 
on FF5. 



FIGURE 11.— Main gear wheel speed during impact on FF5. 


FF5 landing is shown in figure 12 and was derived 
from friction data measurements made by the NASA 
Langley Research Center. 3 

The measured and predicted Orbiter pitch rate on 
FF5 for 1 second after primary gear impact is shown 
in figure 13. The measured pitch rate was sampled at 
a much lower rate than the strain gage data and ap- 
parently had less accurate resolution. Nonetheless, 
the correlation between test and analysis was good. 
The pitch rate transient indicates the initial pitchover 
of the Orbiter after the primary main gear impact. 


FF5 main impact 

horizontal 

velocity 



FIGURE 12.— Main gear tire/ground friction properties during 
spin-up. 



FIGURE 13.— Orbiter pitch rate transient at main gear impact 
on FF5. 


3 Walter B. Horne, “Tire Friction Coefficient Estimates for 
Decelerating and Steering the Space Shuttle Vehicle During Land- 
ing on Conventional, Grooved, and Lakebed Surfaces, NASA 
Langley Research Center Internal Letter, Nov. 5, 1976. 
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Nose gear impact loads on the lower drag brace 
and upper strut trunnion for FF5 are shown in figure 
14, together with the predicted loads using the pre- 
ALT nose gear mathematical model. Only one of the 
Z trunnion loads is shown because the inboard and 
outboard loads are very similar. The correlation be- 
tween analysis and test was fairly good; however, 
refinements in the nose gear unsprung mass pro- 
duced very good correlation as shown in figure 15. 
The remaining differences between analysis and test 
are attributed to the same type inaccuracies as for the 
main gear. The nose gear stroke at nose gear impact 
for FF5 is shown in figure 16. As in the case of the 
main gear, the stroke correlation is excellent during 
initial impact but has an error of about 11 percent at 
the time of maximum stroke. This error is again at- 


tributed to instrumentation, calibration, and model- 
ing inaccuracies. 

The maximum main gear tire loads for each ALT 
flight are shown in table II. Also shown are the cor- 
responding horizontal velocity and elevon position 
for each flight. These loads were calculated from. the 
drag brace and trunnion loads. These tire loads cor- 
respond to a 50-50 distribution of the total ground 
reaction to each tire. There was no instrumentation 
to discern differential tire loads on each gear. The 
predicted tire loads for each flight showed good 
agreement, with a maximum error of 13.8 percent on 
FF5. This type of error was within the tolerance 
band of measurements and flight conditions. As pre- 
viously noted, the tire loads on FF4 were the highest 
of all ALT landings. 




(a) 


200 




(b) (b) 


FIGURE 14. — Nose gear loads and correlation for FF5 using pre- FIGURE 15. — Nose gear loads and correlation for FF5 using up- 

ALT math model, (a) Drag brace, (b) Trunnion (left). dated math model, (a) Drag brace, (b) Trunnion (left). 
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Gear strut stroke s 


Table II . — Maximum Main Gear Tire Loads for ALT 


Flight 



Main gear tire load. N (lb) 



Horizontal 
velocity, 
m/s (knots) 
(a) 

Elevon 

position 

(a) 


Right gear 


Left gear 


Measured 

Predicted 

Measured 

Predicted 

With tail cone on 

1 

311 800 

(70 100 ) 

321 600 

(72 300 ) 

377 600 (75 900 ) 

321 600 

(72 300 ) 

76.1 ( 148 . 0 ) 

- 33 ° 

2 

299 800 

(67 400 ) 

298 000 

(67 000 ) 

— — 

298 000 

(67 000 ) 

70.4 ( 136 . 8 ) 

- 33 ° 

3 

297 600 

(66 900 ) 

325 600 

(73 200 ) 

338 100 (76 000 ) 

325 600 

(73 200 ) 

77.0 ( 149 . 7 ) 

1 

OJ 

o 

With tail cone off 

4 

330 500 

(74 300 ) 

338 100 

(76 000 ) 

366 500 (82 400 ) 

338 100 

(76 000 ) 

83.8 ( 162 . 8 ) 

- 16 ° 

5 

244 200 

(54 900 ) 

278 000 

(62 500 ) 

253 500 (57 000 ) 

278 000 

(62 500 ) 

67.4 ( 131 . 1 ) 

- 12 ° 


a At time of maximum load. 


.42 



CONCLUDING REMARKS 

Correlation between main and nose gear impact 
loads was generally achieved within a tolerance band 
of ±20 percent on load magnitude and ±10 percent 
on response frequency. The differences between 
measured and predicted data are attributed to instru- 
mentation, calibration, and mathematical modeling 
inaccuracies. It should be noted that both magnitude 
and phasing were important in this correlation 


because gear load time histories, rather than just 
maximum load, produce the design-type loading on 
the Orbiter primary structure during main and nose 
gear impact. 

One significant factor in obtaining good landing 
loads was the shunt calibration of the gear strain 
gages after free flight 5. The shunt calibration pro- 
duced significant changes in the gear calibration con- 
stants and thus in the actual gear loads. 

Certain mathematical model updates were made 
because of the approach and landing test load cor- 
relation. These updates included changing the main 
gear unsprung mass from 594.0 to 685.9 kilograms 
(40.7 to 47.0 slugs) and the nose gear unsprung mass 
from 144.5 to 211.6 kilograms (9.9 to 14.5 slugs). The 
unsprung mass as used in the mathematical models 
represents an equivalent mass that reflects accelera- 
tion of the gear and its support structure during spin- 
up/spring-back. Also, the bending damping of the 
main gear strut was increased from 3000 to 6000 
N-s/m (200 to 400 lb-s/ft). 

The general conclusion of this correlation was that 
it provided a good basis for updating and verifying 
main and nose gear loads mathematical models and 
thus provided confidence in predicting Orbiter 
design landing loads. 

Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, September 28, 1979 
986-15-31-05-72 
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APPENDIX 


DERIVATION OF THE EULER TRANSFORMATION 


An Euler transformation transforms vectors from 
a body coordinate system to an inertial coordinate 
system. Consider the origin of the body coordinate 
system (X B , Y B ,Z B ) to be located at the center of mass 
of the vehicle and to be initially alined with the iner- 
tial coordinate system (X h Y h Zj) as shown in figure 
17. (Both coordinate systems are orthogonal.) The Z 7 
axis is directed upward and is parallel to the gravity 
vector, and the A^-F/axes define the plane of the 
landing surface (fig. 17(a)). 

The body coordinate system is oriented with 
respect to the inertial coordinate system by the set of 
Euler angles 9 X , 9 y , and 0 2 . With respect to the vehi- 
cle, 9 X represents the roll angle, 9 V the pitch angle, 
and 9 Z the yaw angle. The Euler transformation will 
be developed for the rotation sequence of roll, pitch, 
and yaw, with each successive rotation taken about 
the indicated body axis. 

The first Euler angle 9 X is developed by rotating 
the body about the X B axis. For this configuration 
(fig. 17(b)), a body prime coordinate system 
(X' b ,Y b ,Z b ) is introduced. The transformation that 
transforms vectors in the body prime system to vec- 
tors in the inertial system is given by 



"l 

0 

0 

lx ' 

0 

cos 6 x 

-sin 6 X 

r 

_0 

sind x 

cos 6 

V’s 


(29) 


(X B ,Y B ,Z B ) is introduced. The transformation that 
transforms vectors in the body double-prime system 
to vectors in the prime system is given by 



0 sin# 

y 

1 0 
0 cos 6 


v " 


(30) 


The third and final Euler angle 9 Z is developed by 
rotating the body about the Z B axis. For this con- 
figuration (fig. 17(d)), a body triple-prime coordinate 
system (X B \Y B \Z B ) is introduced. The transfor- 
mation that transforms vectors in the body triple- 
prime system to vectors in the double-prime system 
is given by 


lx"\ 

\ B 


cos 6 Z 

-sin d z 

0] 

/ v f,, \ 

\ B 1 

Jr" 

y b\ 

L = 

sin d 2 

COS0 Z 

0 

) V 

\ 1 B / 

V'B i 


0 

0 

1 

( 7 nt 1 

\ l b I 


(31) 


Because the triple-prime coordinate system repre- 
sents the final body orientation, then 


The second Euler angle 9 y is developed by rotating 
the body about the F^axis. For this configuration 
(fig. 17(c)), a body double-prime coordinate system 


[ x b\ 


[ x b] 

ytn 

1 B 

= { 

\ Y B 

\ Z Bi 


l z b) 


(32) 
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FIGURE 17. — Euler angle definition, (a) Body and inertial systems coinciding, (b) Rotation about X B axis (0^). (c) Rotation about 
Y' b axis (0^). (d) Rotation about Z' B axis (0p. (e) Final orientation. 
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Combining matrix equations (29) through (32) Matrix equations (30) and (31) yield 

gives 



T BI u T BI u T Bl n 


\ X B 


x b\ 

T BI n T BI 12 T BI 22 


k 

f — 1 

1 1 

II 

y b\ 

J BI 3\ TdI 27 TbI 33_ 


[ Z * 


z bJ 


(33) 


X B = C0S °y { X 'b) + sine yi Z B ) (36) 

X B = cosd z{ X B ) ~ sin9 z{ Y B ) ( 37) 

Y B = + COS0 z( r s) < 38 > 


where T B / U = cos 0 v ,cos 9 Z 

Tb/ 12 =- cos0 y sin0 z 
T B l u = sin 9 y 

Fg/ 21 =sin 0 x .sin 0 V cos 0 z + cos 0 x sin 9 : 
T B , n = — sin 0 x sin 0 v .sin 0 Z + cos 0 x cos 9 : 
t bi 23 =~ sin 9 X cos 9 y 
Tg/ 3l = — cos 0 x sin 0 v cos 0 Z + sin 0 x sin 0 Z 
Tfl /,2 = c °s 0 x sin 0 v sin 0 Z + sin 9 X cos 9 : 
T Bl n = cos 0 x cos 0^, 

where [7^/] is the Euler transformation that 
transforms vectors expressed in the body system to 
vectors in the inertial system. 

The Euler angle transformation [T BI ] is 
orthogonal and thus its inverse is equal to its 
transpose. 


Z B Z B 

Substituting equations (36) through (39) into equa- 
tion (35) gives 

^x{^b) + ^(^ 5 ) + ( Z b) 

= (d x cos d y cos 6 z + d y sin d^ (x^ 

+ (o y cos e z - e x Sin d z cos djj (y^J 

+ (e x sin d y + e z ) (z B ) (40) 


M" 1 = I t biT - I t ib] < 34 > 


Putting equation (40) in matrix form 


During a dynamic situation, the Euler angles are 
changing with time such that the integral of the Euler 
rate vector yields the Euler angles. The Euler rate 
vector is the projection of the body rate vector (co x , 
co p co z ) (expressed in the body coordinate system) on 
the nonorthogonal Euler rate axes system 
(X'b'Y'b'Zb) or, mathematically, 

[«A) + «,(>*) + «, ( Z s) 
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1 X 1 

“y 

1 1 

II 
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w z 
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cos 0 cos 0 sin 0 0 
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y z z 



-sin 0 z cos d y cos 6 Z 0 



_ sin dy 0 1_ 




( 41 ) 
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Solving matrix equation (41) for the Euler rate 
vector (ff x , 9 y , 0 Z ) 




CO ) 

z J 


or 


0 = 
X 




co x cos 6 z 


u y sin e z 


26 2 ° y + e y“z) 


cos 0„ 


(44) 


= sin 6 Z + co cos 6 Z + 6 x 6 z cos 6 (45) 


d z = -6 X sin Q y - e y 6 x cos 6 y + co z 


(46) 


cos 0 Z 

—sin 0 z 

0“ 

K 

sin 6 Z 

cos 0 Z 

0 

Ico 

X y 

-sin 0„ cos 0 

sin 0 sin 0^ 

1 

JO) 

L y z 

y z 


\ z l 


(42) 


where the matrix [T BE ] is the inverse of [T ES \. 

The time derivative of equation (42) yields an ex- 
pression for the time rate of change of the Euler rate 
components. 


Differential equations (44) through (46) express the 
body angular rate and angular acceleration compo- 
nents in a form that can be conveniently integrated 
to yield the angular coordinates 0 X , 0 G z (the three 
Euler angles) that define the orientation of the rigid 
vehicle with respect to the inertial coordinate 
system. Many numerical integration schemes are 
available for solving equations (44) through (46). 
The numerical integration scheme used in the land- 
ing simulation was a variable step, with the time step 
set by the truncation error. This type of integration 
scheme worked satisfactorily for the Shuttle landing 
dynamics simulation. 
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